Variational description of multi-fluid hydrodynamics: Uncharged fluids. 
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We present a formalism for Newtonian multi-fluid hydrodynamics derived from an unconstrained 
variational principle. This approach provides a natural way of obtaining the general equations of 
motion for a wide range of hydrodynamic systems containing an arbitrary number of interacting flu- 
ids and superfluids. In addition to spatial variations we use "time shifts" in the variational principle, 
which allows us to describe dissipative processes with entropy creation, such as chemical reactions, 
friction or the effects of external non-conservative forces. The resulting framework incorporates 
the generalization of the entrainment effect originally discussed in the case of the mixture of two 
superfluids by Andreev and Bashkin. In addition to the conservation of energy and momentum, 
we derive the generalized conservation laws of vorticity and helicity, and the special case of Ertel's 
theorem for the single perfect fluid. 

We explicitly discuss the application of this framework to thermally conducting fluids, superfluids, 
and superfluid neutron star matter. The equations governing thermally conducting fluids are found 
to be more general than the standard description, as the effect of entrainment usually seems to 
be overlooked in this context. In the case of superfluid 4 He we recover the Landau-Khalatnikov 
equations of the two-fluid model via a translation to the "orthodox" framework of superfluidity, 
which is based on a rather awkward choice of variables. Our two-fluid model for superfluid neutron 
star matter allows for dissipation via mutual friction and also "transfusion" via /9-reactions between 
the neutron fluid and the proton-electron fluid. 



I. INTRODUCTION 

The main purpose of this work is to develop a formal- 
ism that allows one to derive the equations of motion 
for a general class of multi-constituent systems of inter- 
acting charged and uncharged fluids, such as conducting 
and non-conducting fluids, multi-fluid plasmas, superflu- 
ids and superconductors. For the sake of clarity of pre- 
sentation we restrict ourselves here to uncharged fluids, 
while the case of charged fluids and their coupling to the 
electromagnetic field will be treated in a subsequent pa- 
per [1]. 

Long after the completion of classical Hamiltonian par- 
ticle mechanics, the quest of finding a variational (or 
"Hamlitonian") description of hydrodynamics has surpris- 
ingly been a long-standing problem, which started only 
a few decades ago to be fully understood. The reason 
for this can be traced to the nature of the hydrodynamic 
equations, which are most commonly expressed in their 
Eulerian form in terms of the density p and velocity v, 
where the information about the underlying flowlines has 
been hidden. Fluid particle trajectories, i.e. flowlines, 
can still be recovered by integrating the velocity field, 
but they are not independent quantities of the Eulerian 
description. However, it turns out that the "true" fun- 
damental field variables of Hamiltonian hydrodynamics 
are the flowlines, which determine p and v as derived 
quantities. 

Consider as an example the Lagrangian density A de- 
scribing a barotropic perfect fluid, which in analogy to 



classical mechanics one would postulate to be 



A(p,v) 



^ P v 2 -£(p), 
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where £(p) represents the internal energy density of the 
fluid. We note that the internal energy defines the chem- 
ical potential p and the pressure P as 

d£ = p dp , and P + £ = pp . 

The corresponding action is defined in the usual way as 
1= J AdVdt, and the variation SA of the Lagrangian 
density is 

5A = pv ■ Sv + (v 2 /2 -p)5p. 

Requiring the action X to be stationary with respect to 
free variations Sp and Sv is immediately seen to be use- 
less, as this leads to the over-constrained equations of mo- 
tion pv = and p — v 2 /2. In fact, it has been shown [2] 
that an unconstrained variational principle with p and v 
as the fundamental variables cannot produce the Eulerian 
hydrodynamic equations. The reason for this is rather in- 
tuitive, as it is evident that free variations of density and 
velocity probe configurations with different masses (i.e. 
different numbers of particles) , which is not an actual de- 
gree of freedom of the dynamics of the system. Therefore 
the variational principle has to be constrained or refor- 
mulated in some way in order to restrict the variations 
to the physically meaningful degrees of freedom. 

The historic approach to this problem in Newtonian 
physics has been to supplement the Lagrangian with ap- 
propriate constraints using Lagrange multipliers. This 
method was pioneered by Zilsel [3] in the context of the 
two-fluid model for superfluid 4 He, who used the con- 
straints of conserved particles (i.e. mass) and entropy. 
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However, as pointed out by Lin [4], this is generally in- 
sufficient, as it results in equations of motion restricted 
to irrotational flow in the case of uniform entropy. Lin 
showed that one has to add a further constraint, namely 
the "conservation of identity" of fluid particles in order 
to obtain the most general hydrodynamic equations. We 
can label particles by their initial positions a, and so we 
can write their flowlines (a, t). The famous "Lin 

constraint" is dta + v ■ Va = 0, i.e. the identity or label 
of a particle is conserved under its transport. For reviews 
of this approach and its relation to the "Clebsch repre- 
sentation" we refer the reader to [5-7], and references 
therein. 

Although this method produces the correct equations 
of motion, it does not seem very natural due to the rather 
ad hoc introduction of constraints, and the need for un- 
physical auxiliary fields (the Lagrange multipliers). It 
was pointed out by Herivel [8] that the Lagrangian as op- 
posed to Eulerian formulation of hydrodynamics results 
in a much more natural variational description, and this 
approach was further developed and clarified by Seliger 
and Whitham [5]. Instead of using p and v as funda- 
mental variables, hydrodynamics can also be understood 
as a field theory in terms of the flowlines x(a,t), or 
equivalently a = a(x,t). It turns out that this formu- 
lation allows for a perfectly natural unconstrained varia- 
tional principle. This seems rather intuitive considering 
that hydrodynamics is a smooth-averaged description of 
a many-particle system, which is described by a varia- 
tional principle based on the particle trajectories, i.e. xn 
and xpf. 

We can express the velocity and density in terms of the 
flowlines as v — dtx(a,t) and p(x,t) = po(a)/ det(J r %), 
where J l j — dx l jdo? is the Jacobian matrix correspond- 
ing to the map a i— > x(a, t) between the physical space x 
and the "material space" a. Any further comoving quan- 
tities like the entropy s are determined in terms of their 
initial value so (a). Substituting these expressions into 
the Lagrangian A, one obtains an unconstrained varia- 
tional principle for the field x(a,t), which results in the 
correct equations of motion. It is interesting to note that 
this approach implicitly satisfies Lin's constraint, as we 
are varying the particle trajectories x(a,t), along which 
a is a constant by construction. Also, we do not need 
to impose an a priori constraint on the conservation of 
mass, as it is automatically satisfied by these "convec- 
tive" variations: shifting around flowlines obviously con- 
serves the number of flowlines, and therefore the number 
of particles. One can actually derive the Lin constraint 
by transforming this Lagrangian framework back into a 
purely Eulerian variational principle [5, 6], which shows 
that these two approaches are formally equivalent. 

As pointed out by Bretherton [9], one can even more 
conveniently use a "hybrid" approach, in which the La- 
grangian is expressed in terms of the Eulerian hydro- 
dynamic quantities v, p, s etc, but one consider them 
as functions of the underlying flowlines. Their varia- 
tions are therefore naturally induced by variations £ of 



the flowlines x(a,t). In general relativity the same idea 
was pioneered by Taub [10], and has subsequently been 
largely developed and extended by Carter [11-13], who 
also coined the term "convective variational principle" for 
this approach. Carter and Khalatnikov [14] have further 
demonstrated the formal equivalence of the convective 
approach and the more common Clebsh formulation that 
results from an Eulerian variational approach. A "trans- 
lation" of the covariant convective formalism into a New- 
tonian framework (albeit using a spacetime-covariant lan- 
guage close to general relativity) is also available [15, 16]. 
The convective approach in relativity has independently 
been developed by Kijowski [17], and Hamiltonian for- 
mulations have been constructed by Comer and Langlois 
[18] and Brown [19]. Here we are using the convective 
(or "hybrid") variational principle in order to derive the 
Newtonian multi-fluid equations, and our notation and 
formalism follows most closely the framework developed 
by Carter. 

We conclude our example of the simple barotropic fluid 
by using the convective variational principle to derive the 
Euler equation. The expressions for (Eulerian) variations 
of density and velocity induced by infinitesimal spatial 
displacements £ of the flowlines are well known 1 (e.g. see 
[20]), namely 

Sp = -V (p£) , and Sv = dt£ + (v • V)£ - (£ • V)t» . 

Inserting these expressions into the variation of the ac- 
tion 61 = J SA dV dt with 5 A given above, and after some 
integrations by parts and dropping total divergences and 
time derivatives (which vanish due to the boundary con- 
ditions), we find 
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p(d t + v ■ V)v + pVp 



+v {dtp + V • (pv)} 



dVdt. 



If we assume conservation of mass 2 , i.e. 
dtp + V • (pv) = 0, then stationarity of the action 
(i.e. SI = 0) under free variations £ directly leads to 
Euler's equation, namely 

(d t + v ■ V) v + -VP = , 
P 

where we have used the thermodynamic identity 
pVp, = VP. This shows that an unconstrained convec- 
tive variational principle produces to the correct hydro- 
dynamic equations of motion in a surprisingly simple and 
straightforward way. 



1 A generalization of these expressions to include time-shifts is 
derived in Appendix A 

2 This will be seen to be a consequence of the variational principle 
rather than an a-priori assumption when time-shift variations are 
included. 
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The spatial variations £ have three degrees of freedom, 
resulting in one vector equation, which represents the 
conservation of momentum. In order to complete the 
description we will need a fourth variational degree of 
freedom to produce the missing energy equation. This 
can be achieved by considering time-shifts, which are a 
natural part of the covariant relativistic approach, but 
which we have to be considered explicitly in the conven- 
tional "3+1" language of Newtonian space-time. These 
time-shifts variations allow us to take this formalism to 
its full generality, as we can now describe even dissipative 
processes with entropy creation, particle transformations 
(i.e. chemical reactions), resistive frictional forces etc. 
These dissipative systems are of course still conservative 
as long as one includes entropy, which is why they can be 
described by an action principle. The second law of ther- 
modynamics, however, is obviously not contained in the 
action principle and has to be imposed as an additional 
equation on the level of the equations of motion. 

We note that the equations we derive here do not ex- 
plicitly include shear- and bulk-viscosity effects. How- 
ever, the current form of the equations is in principle 
general enough to allow for both of these effects: bulk 
viscosity is caused by heat flow or chemical reactions due 
to thermal or chemical disequilibrium, both of which can 
already be described in the current formulation. Shear 
viscosity on the other hand has to be introduced as an 
"external" force, the problem therefore consists in pre- 
scribing a physically reasonable model for a multi-fluid 
generalization of the shear stresses. Including viscosity 
should therefore not be a matter of actually extending 
the current framework but rather of appropriately ap- 
plying it in order to describe such processes. An explicit 
discussion of this is postponed to future work. Further 
work is also necessary in order to extend this formal- 
ism to include elasticity (as pioneered in the relativis- 
tic framework [21]), and especially to allow for an elas- 
tic medium interpenetrated by fluids as encountered in 
the inner neutron star crust, or any type of conducting 
solid. As shown in [22], a Kalb-Ramond type extension 
is required for the macroscopic treatment of quantized 
vortices in superfluids. With the present formalism we 
can describe superfluids either on the local irrotational 
level, or on the smooth-averaged macroscopic level by 
neglecting the (generally small) anisotropy induced by 
the quantized vortices. 

The plan of this paper is as follows: in Sect. II we de- 
rive the general form of the equations of motion for multi- 
constituent systems using the convective variational prin- 
ciple. In Sect. Ill we show the conservation of energy 
and momentum implied by these equations. In Sect. IV 
we derive conserved quantities under transport by the 
flow, namely the vorticity and helicity. We then give 
the explicit functional form of the Lagrangian density 
for hydrodynamic systems in Sect. V, and in Sect. VI we 
discuss several applications of the foregoing formalism to 
particular physical systems. 



II. VARIATIONAL DESCRIPTION OF 
MULTI-CONSTITUENT SYSTEMS 

A. Kinematics 

We want to describe systems consisting of several con- 
stituents distinguished by suitably chosen labels, and we 
use capital letters X, Y , ... as indices which run over these 
constituents labels. As the fundamental quantities of the 
kinematic description we choose the constituent densities 
nx and the associated transport currents nx, which are 
related to the respective velocities vx as 

nx = nxvx , where X e {constituent labels} . (1) 

Not all constituents can necessarily move independently 
from each other, i.e. not all velocities Vx have to be 
different: viscosity and friction due to particle collisions 
on the microscopic level can effectively bind constituents 
together on very short timescales. We therefore distin- 
guish between the notions of constituents X , characteriz- 
ing classes of microscopic particles, and fluids, which are 
sets of constituents with a common velocity. 

We note that in this framework entropy can be de- 
scribed very naturally as a constituent for which we re- 
serve the label X = s, and we write 



(2) 



where s is the entropy density. In this context it is in- 
structive to think of the entropy as a gas of particle-like 
thermal excitations (e.g. phonons, rotons etc.), which 
makes its description as a constituent on the same foot- 
ing with particle number densities quite intuitive. 



B. Dynamics 

The dynamics of the system is governed by an action 
J defined as 



/ 



A H dV dt , 



(3) 



in terms of the hydrodynamic Lagrangian A H . The La- 
grangian density A H depends on the kinematic variables, 
which are the densities nx and the currents nx, i.e. 
A H = A H (nx, nx ). The total differential of A H defines 
the dynamical quantities p§ ("energy") and jf 1 ("momen- 
tum") per fluid particle as the canonically conjugate vari- 
ables to nx and nx, namely 

dA u =^2(p^dn x +p < - dn x ), so ^o = ^~' p c =^-, 

"(4) 

where here and in the following the sum over repeated 
constituent indices is explicitly indicated by a E, i.e. no 
automatic summation convention applies to constituent 
indices. 
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C. The convective variational principle 

As we have seen in the introduction, one cannot ap- 
ply the standard variational principle to A H in terms of 
the Eulerian hydrodynamics variables nx and % . From 
(4) it is obvious that allowing free variations of densities 
Snx and currents Snx would lead to the trivial equa- 
tions of motion p§ — and = 0. Instead, we con- 
sider the Lagrangian to be a functional of the underly- 
ing flowlines ar^ = ar^ (a* , i) , and therefore admit only 
variations 8n x , dn x that are induced by infinitesimal dis- 
placements of the flowlines. These "convective" variations 
naturally conserve the number of particles (i.e. the num- 
ber of flowlines) and no constraints are required in the 
variational principle as was discussed in more detail in 
the introduction. 

We apply infinitesimal spatial displacements £ x and 
time-shifts Tx to the flowlines of the constituent X. The 
resulting induced variations of density and current have 
been derived in Appendix A, namely the density varia- 
tion (A18) for constituent X is 

5n x = -V • [nx€ x ] + \ n x ■ Vt^ - T X d t n x ] , (5) 

while the current variation Snx is given by (A20) and 
reads as 

8n x = n x dt£x + ( n x ■ V) £ x ~ i€x ' v ) n * 

-nx(V-tx)-dt(nxVc) ■ (6) 

Inserting these expressions into the variation of the La- 
grangian (4) and integrating by parts, we can rewrite the 
induced variation 5A H in the form 

SA H =J2(9 X r x - f ■ix)+d t R + ^ R- (7) 

The time derivative and divergence terms will vanish in 
the action integration (3) by the appropriate boundary 
conditions (i.e. £ = and r = 0) and are irrelevant 
as far as the variational principle is concerned, but for 
completeness we note that their explicit expressions are 

R = ( nx ^ P* 7 *) - ( 8 ) 

-«x(%i4+%y)]. (9) 

The induced action variation therefore has the form 

51 = J2j(9 X T X -f X - Sx) dVdt , (10) 

where the force densities f x (acting on the constituent) 
and the energy transfer rates g x {into the constituent) 
are found explicitly as 

f X = n x {d t ^ -X7p*)-n x x (Vx/)+//x, (11) 

g x = vx-{f -^r x )-^r x , (12) 



where Fx is the particle creation rate for the constituent 
X, i.e. 

r x = d t n x + V • n x . (13) 

The force density f x is the total momentum change rate 
of the constituent X, and we see that the last term in 
(11), i.e. the "rocket term" j^Tx, represents a contri- 
bution that is purely due to the change of the particle 
number. Therefore it will be convenient to define the 
purely "hydrodynamic force" f% , as 

^En^V-V^-nxxlVx}/). (14) 

With this definition we can now write the force density 
(11) and energy transfer rate (12) in the form 

f X = f^+^fx, (15) 
/ = vx-fi-^rx. (16) 

D. The equations of motion 

Up to this point we have developed only purely math- 
ematical identities without a specific physical content. 
The equations of motion are obtained by imposing which 
type of invariance the action X should satisfy under cer- 
tain infinitesimal variations. The most general equations 
are obtained by requiring that a common displacement 
£x = £ an d time shift Tx = t of all constituents should 
result in an action variation of the form 

51= J (<?ext r - / cxt • dV dt , (17) 

where / ex t and g ex t are interpretable as the external force 
density and energy transfer rate. This generalizes the 
more common action principle of isolated systems, in 
which the external influences / cx t and g cx t vanish and 
therefore the equations of motion are obtained by requir- 
ing the action to be invariant under small variations. 
"External" here is meant in the sense of not being in- 
cluded in the total Lagrangian, which could also mean, 
for example viscous or gravitational forces. The resulting 
minimal equations of motion obtained from comparing 
with (10) are therefore found as 

^/ X = /cxt, and 2$*= Se*t- (18) 

Together with (11) and (12) this represents the Euler- 
Lagrange equations associated with this variational prin- 
ciple. If all constituents X form a single fluid, namely all 
constituents have a common velocity, then only common 
displacements of all constituents make sense in the vari- 
ational principle. For this class of non-conducting mod- 
els, (18) represent the full equations of motion obtainable 
from the variational principle. In order to complete the 
model, one has to specify the hydrodynamic Lagrangian 
A H , the external interactions / ex t and g ext , and the cre- 
ation rates Fx as functions of the kinematic variables. 
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In the case of conducting models, at least some of the 
constituents are allowed to move independently, the sys- 
tem therefore consists of more than one fluid. This in- 
creases correspondingly the number of degrees of free- 
dom, and more equations of motion are required. They 
are obtained very naturally from the variational prin- 
ciple, as independent displacements (in space and time) 
are permitted for each fluid. Therefore the resulting force 
acting on each fluid can be prescribed by the model, sub- 
ject to the restriction only of satisfying the minimal equa- 
tions of motion (18). 

As an example, consider the case of a simple conduct- 
ing model consisting of two fluids, where we useX and Y 
are constituent indices running only over the respective 
constituent labels, i.e. X S {fluid 1} and Y € {fluid 2}. 
We then have the respective force densities acting on each 
of the two fluids as = J2x f X and /(2) = EJ Y . 
which by (18) have to satisfy /(i) + f(2) = /ext- There- 
fore there are now exactly two force densities (e.g. fm 
and /ext) freely specifiable in the model, corresponding 
to the additional degrees of freedom of two fluids. In this 
case /(i) could for example represent a mutual force the 
two fluids exert on each other, e.g. a resistive friction 
force. 



III. CONSERVATION OF ENERGY AND 
MOMENTUM 

Using the explicit expression (11) for the force density 
f x , we can write 

-^(nxVp*+4Vp*) . (19) 

We define the "generalized pressure" \P via the total Leg- 
endre transformation of A H , namely 

* = A H - ^ (n x p£ + -p*) , (20) 

which is seen from (4) to result in the total differential 

rf* = - ( n x d P% + n x ■ dp*) > (21) 

and therefore the last sum in (19) is simply V^. We 
can now cast the force equation (18) in the form of a 
conservation law for the total momentum, namely 

dtJt + V^ =f ext , (22) 

where the hydrodynamic momentum density J H and 
stress tensor are defined as 

J H = n xP" , and TS = £ r&p* ] + * ^ , (23) 

and where gij are the components of the metric tensor de- 
termining the relation between physical distance dl and 



coordinate intervals dx % , i.e. dl 2 = gij dx % dxJ . In Carte- 
sian coordinates this is simply — 5ij. A proof of 
the symmetry of the stress tensor together with a 
more elegant derivation of momentum conservation as a 
Noether identity of the variational principle is given in 
Appendix B. 

Using expressions (11) and (12), we can further show 
that 

J2^ X = J2 i nx ■ dt ^ - nx ■ v ^ - rx $] 

-^(pSdtnx+T? -dtnx), (24) 

and we see from (4) that the last sum simply represents 
<9tA H . We can therefore rewrite the energy equation (18) 
in the form of a conservation law, namely 

d t E a + V • Q H = 5cxt , (25) 

where the hydrodynamic energy density E H and energy 
flux Q H are given by 

Ev = J2 n x-^ - A h> and g H = ^(-pJ)n x . (26) 

We see that the energy density E n has quite naturally 
the form of a hamiltonian, i.e. H H (nx,p x ) = E H , as it is 
the Legendre-transformed (with respect to the momenta) 
of the Lagrangian A H . 

IV. CONSERVATION ALONG FLOWLINES 

In addition to the total energy-momentum conserva- 
tion, derived in the previous section, we can find fur- 
ther conserved quantities for individual constituents, for 
which conservation holds under transport by the fluid 
flow. Because the following derivations apply to individ- 
ual constituents instead of the sum over all constituents, 
we will omit the constituent index X in this section in 
order to simplify the notation. 

Transport of a quantity by the fluid flow is closely re- 
lated to the Lie derivative with respect to the fluid ve- 
locity, therefore these conservation laws are most easily 
derived using the language and theorems of differential 
forms instead of vectors. We will use this formalism in 
deriving the transport-conservation laws, but we also give 
the essential steps and results translated in the more com- 
mon vector- and index-notation, so that familiarity with 
exterior calculus should not be necessary (albeit helpful) 
for reading this section. 

A. Kelvin-Helmholtz vorticity conservation 

We define the vorticity 2-form w (with components 
Wij ) as the exterior derivative (denoted by d) of the mo- 
mentum 1-form p (with components Pi), namely 

w = dp, i.e. w i3 = 2V[iPj] , (27) 
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where [ij] denotes antisymmetric averaging, i.e. 
2A[iBj] — AiBj — AjBi. In three dimensions we can de- 
fine the more common vorticity vector W as the dual 
(with respect to the volume form tijk) of the 2-form w, 
namely 

VP ee \e iik w jk = (V x p) 1 . (28) 
The volume form is defined as 

tijh = Vg[i,j,k] , (29) 

where g = det(gij) and k] is the sign of the permu- 
tation of {1, 2, 3}, which is zero if two indices are equal. 
The duality between w and W implies 

Wij = e ijk W k , (30) 

which is easily verified by inserting (28). We note that 
due to the Poincare property (namely dd = 0) , the exte- 
rior derivative of the vorticity 2-form vanishes identically, 
i.e. 

dw = V W = . (31) 

We can rewrite the expression (14) for the hydrodynamic 
force / H in the language of forms as 

d t p + v\dp-dp = ^-l ii , (32) 

where J indicates summation over adjacent vector- and 
form- indices, i.e. in this case (v\dp)i — 2v 3 VuP{\. In the 
following it will be convenient to separate the force per 
particle into its non-conservative part J an d a conserva- 
tive contribution d<p, namely 

-/ H = d0 + $. (33) 

1 1 > Cartan formula for the Lie derivative of a p-form 
applied to the 1-form p yields 

£ v p = v\dp + d(v\p), (34) 

which in explicit index notation reads as 
£ v pi = 2v 3 VyPi] + Vi(^pj). Using this identity 
and (33) we rewrite the force equation (32) more 
conveniently as 

(fit + £ v )p = dQ + 3, (35) 

where the scalar Q is given by Q = pa + v\p. Lie deriva- 
tives and partial time derivatives commute with exterior 
derivatives, so we can apply an exterior derivative to (35) 
and obtain the Helmholtz equation of vorticity transport, 
namely 

(dt + £ v )w = dS, (36) 

which shows that the vorticity is conserved under trans- 
port by the fluid if and only if the hydrodynamic force 



per particle acting on the fluid is purely conservative, i.e. 
if 5 = 0. In its more common dual form, this equation 
can be written as 

d t W -Vx(»xf)=Vx5. (37) 

The Helmholtz vorticity conservation expresses the 
conservation of angular momentum of fluid particles, 
and we can equivalently derive it in its integrated form, 
namely the conservation of circulation as first shown by 
Kelvin. We consider a 2-surface E and define the circu- 
lation C around its boundary fiE as 

C = <f> p=f Ptdx i . (38) 

Using Stoke's theorem, we see that the circulation around 
fiE is equivalent to the vorticity flux through the surface 
S, i.e. 

C = [ w = \ ( Wii dx l A dx j , (39) 

and the more familiar dual expression is found by insert- 
ing (30): 

C = J W-dS, (40) 

where the surface normal element dS is 
dSi = \eijkdx 3 l\dx k . Using (35) the comoving 
time derivative of the circulation C yields 

§ = i/«*-/<* + (41) 

which is known as Kelvin's theorem of conservation of 
circulation. As we have already seen before, strict con- 
servation only applies if the non-conservative force per 
particle J vanishes. 



B. Vorticity and superfluids 

The hydrodynamics of superfluids is characterized by 
two fundamental properties: on one hand by the absence 
of dissipative mechanisms like friction or viscosity, and on 
the other hand by irrotational flow. As we will see now, 
the hydrodynamic description of superfluids is therefore 
a natural subclass within the more general framework 
of multi-constituent hydrodynamics presented here. Let 
us assume that a constituent X — S is superfluid, with 
particle density ns, velocity i>s and mass m s . The ab- 
sence of microscopic dissipative mechanisms implies that 
the superfluid is not bound to any other constituents , 
i.e. it is a perfect conductor in the sense that it can 
flow freely even in the presence of other constituents. 
Dissipation-free flow is characterized by the absence of 
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non-conservative forces acting on the bulk 3 of superfluid, 
i.e. 

3 s = . (42) 

As a consequence of (36) and (41) we see that the vortic- 
ity (and therefore circulation) of a superfluid is strictly 
conserved. The second constraint, which distinguishes 
a superfluid from a perfect fluid, is that a superfluid is 
locally irrotational, i.e. its vorticity is zero, so 

w s = , <^=> W s = . (43) 

Due to the vorticity conservation of superfluids, this con- 
straint remains automatically satisfied if it is true at some 
instant t, i.e. it is consistent with the hydrodynamic evo- 
lution. 

The formulation most commonly found in the litera- 
ture on superfluids and superconductors is based on the 
concept of the so-called "superfluid velocity", which is 
constrained to be irrotational [23, 24]. If one interpreted 
this as the actual transport-velocity vg, such a constraint 
would generally not be consistent with the equations of 
motion, contrary to the natural conservation of the mo- 
mentum vorticity w S • This "orthodox" formulation of su- 
perfluidity, which goes back to Landau's two-fluid model 
for 4 He, is therefore a rather unfortunate misinterpreta- 
tion of physical quantities, as the so-called "superfluid 
velocity" is necessarily to be interpreted as the rescaled 
superfluid momentum in order to make this constraint 
consistent with hydrodynamics. The fact that in New- 
tonian single-fluid contexts the particle momentum only 
differs by a constant mass factor from the velocity has un- 
fortunately lead to a less than careful distinction between 
these fundamentally different quantities. This simple 
identification no longer holds true in more general con- 
texts, like in the case of multi-fluids (e.g. superfluids) or 
even in the case of a single relativistic perfect fluid. The 
velocity-circulation is generally not conserved, contrary 
to the conservation of momentum circulation (41). The 
orthodox framework of superfluid hydrodynamics will be 
discussed in more detail in Sect. VI D. 

In addition to the superfluid constraints of being 
dissipation-free and irrotational, there is a further impor- 
tant restriction, namely the quantization of circulation. 
An irrotational flow can still carry non-zero circulation 
in the presence of topological defects (such as vortices). 
In order to see this, we note that (as a consequence of 
(43)) we can write the superfluid momentum p as the 
gradient of a phase tp, namely 

p s = hdip, i.e. p s = h\7tp. (44) 



3 However, there can be a non-conservative force acting on the 
superfluid at a vortex-core if the vortex is pushed by another 
fluid. This mechanism gives rise to the so-called effect of "mutual 
friction". 



The circulation (38) can therefore be non-zero if en- 
closes a topological defect in ip, i.e. a region where ip 
(and p s ) is not defined, as for example in the case of flow 
inside a torus. While in the case of a perfect irrotational 
fluid the resulting circulation could have any value, the 
superfluid phase <p is restricted to change only by a mul- 
tiple of 2ir after a complete tour around the defect. The 
resulting circulation is therefore quantized as 

C = 2Nnh, with N £ Z , (45) 

which gives rise to the well-known quantized vortex struc- 
ture of superfluids. 

C. Helicity conservation 

Contrary to the conservation laws derived in the pre- 
vious sections, which have been known for more than a 
century, there is a further conserved quantity namely the 
so-called helicity, whose existence in hydrodynamics has 
only been pointed out comparatively recently by Moffat 
[25]. This quantity is analogous to the magnetic helicity 
conservation found in magneto-hydrodynamics [26], and 
it is related to the topological structure of the vorticity, 
i.e. its "knottedness" [27]. The relativistic analogue of 
this conservation has been shown by Carter [13, 28, 29], 
and generalizations have been discussed by Bekenstein 
[30]. 

We define the helicity 3-form H_ (with components 
Hijk) as the exterior product of the momentum 1-formp 
with the vorticity 2-form w, i.e. 

H=pAw, (46) 

which in components reads as H%jh — 3p[iWjk] • A 3-form 
in a 3-dimensional manifold is dual to a scalar, so we can 
define the helicity density h as 

Hijk = h Cijk ■ (47) 

From the duality relation together with the definition 
(46), we see that the helicity scalar has the following 
explicit expression 

h=j^ k H ijk = Pi ^ k w jk =p\W=p-(Vxp). (48) 

Using (35) and (36), the comoving time-derivative of H 
can be expressed as 

(d t + £ v )H = [{d t + £ v )p\ Aw+pA[{d t + £ v )w\ 

= (dQ + $) Ai + pAd| 

= d(Qw) + [d(p A I) + 2dl A p] . (49) 

We see that, not surprisingly, the vanishing of the non- 
conservative force J is a necessary (albeit not sufficient) 
condition for the conservation of helicity. We introduce 
the total helicity H of a volume V as 

H= f H= f hdV , (50) 
Jv Jv 
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and, assuming 5 = 0, we find for the comoving time 
derivative of H: 



Qw - 



av 



av 



QW-dS. (51) 



The helicity Ti. of a volume V is therefore conserved un- 
der transport by the fluid if, in addition to 5 = 0, the 
vorticity W vanishes on the surface dV surrounding this 
volume. 



V. HYDRODYNAMICS 

A. The Lagrangian of hydrodynamics 

In the previous sections we have derived the most gen- 
eral form of the Euler-Lagrange equations (18) associ- 
ated with the convective variational principle, together 
with the force densities (11) and energy transfer rates 
(12). We are now interested in a particular class of La- 
grangian densities A H , namely those which describe New- 
tonian hydrodynamics. One can postulate the general 
form of the hydrodynamic Lagrangian A H in analogy to 
canonical particle mechanics as 



A H {n x ,n x ) = Vm x ^ -£ 



(52) 



where £ is a thermodynamic potential related to the in- 
ternal energy (or "equation of state") of the system. We 
therefore find the following general form for the conjugate 
momenta p% and as defined in Eq. (4) : 



2 on x 



d£_ 

dn- 



■x 



(53) 



We want to identify these conjugate momenta with the 
actual physical energy and momentum per fluid particle, 
which implies that under a Galilean boost —V inducing 
the transformations 

v' x = vx + V , n' x =n x , d t ' = d t - V ■ V , (54) 

these momenta should transform (e.g. see [23, 31]) as 

' = -J^ + V-t^ + ^m x V 2 , and p*' =p K +m x V. 

(55) 

One can verify that in this case the hydrodynamic force 
densities f% defined in (14) are invariant under Galilean 
boosts as one should expect. The particle creation rates 
Fx defined in (13) are also Galilean invariant, so that the 
transformation of the total force densities f x of (11) is 
seen to be 



f x ' = f x + Vm x r x . 



(56) 



The equations of motions of an isolated system, i.e. 
^2 f x = 0, are therefore Galilean invariant if and only 
if the total mass is conserved, i.e. if 



^m x r x =0. 



(57) 



By using (55) we can show that the energy transfer rates 
(16) transform as 



V 2 



x ■ 



2 ' 



(58) 



and due to mass conservation (57) the total energy 
change rate therefore satisfies 



(59) 



so that the total energy conservation of an isolated sys- 
tem is Galilean invariant. 

In general the transformation properties (55) are only 
consistent with the conjugate momenta (53) if £ is itself 
Galilean invariant, which is shown in Appendix C. This 
implies that the velocity dependence of £ can only be of 
the form 



£{n x ,n x ) = £(ux,Axy) , 



(60) 



where Ax y is the relative velocity between fluid X and 
fluid Y, i.e. 



A x Y = v x - v Y = 



Tlx fly 



Tlx 



71 y 



(61) 



We note that a function £ of the form (60) satisfies the 
identity 



x . d£ 



0. 



(62) 



which can be used together with (53) to show that the 
hydrodynamic momentum density (23) satisfies 



J H = E = E mXn x = P i 



(63) 



i.e. the hydrodynamic momentum density «7 H is equal 
to the total mass current p as a consequence of Galilean 
invariance. 

In addition to the requirement of Galilean invariance 
we will restrict our attention to systems of "perfect" 
multi-constituent fluids in the sense that their energy 
function £ is isotropic. This means that we consider only 
equations of state of the form 



£(n x ,A XY ) = £(n x ,A XY ). 



(64) 



Summarizing we can now write the hydrodynamic La- 
grangian density (52) for this class of perfect multi-fluid 
systems as 

A n (n x ,n x ) = y2 m x ^-£(n x ,A XY ). (65) 
^ 2n x 

It is interesting to note that contrary to the relativis- 
tic case, which is governed by a fully covariant hydrody- 
namic Lagrangian density (e.g. see [13]), the Newtonian 
Lagrangian (65) is not strictly Galilean invariant because 
of the kinetic energy term. The violation is sufficiently 
weak, however, that is does not affect the Galilean in- 
variance of the resulting equations of motion. 
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B. Conjugate momenta and entrainment effect 

The total differential of the energy function 
£(n x , A XY ) represents the first law of thermodynamics 
for the given system, namely 

d£ = ^^ dn x + \^2<x XY dA XY , (66) 

X,Y 

which defines the chemical potentials ^ and the sym- 
metric entrainment matrix a* Y as the thermodynami- 
cal conjugates to nx and A XY . The conjugate momenta 
(53) are therefore explicitly found as 

P* = rr^vx-Y] A XY , (67) 

-p* = ^-m x ^-+v x -p*- (68) 

The expression (67) for the momenta in terms of the ve- 
locities is interesting, as it shows that in general the mo- 
menta are not aligned with the respective fluid velocities, 
which is the so-called entrainment effect 4 . The simple 
single-fluid case, in which the momentum is just p = mv, 
is only recovered if there is no entrainment between the 
fluids (i.e. a XY = 0) or if all constituents move together 
(i.e. Axy = 0). This phenomenon is well-known (albeit 
not under the name entrainment) in solid-state physics, 
for example the electron momentum in a crystal lattice 
is connected to its velocity by an effective mass-tensor 
(e.g. see [33]). For a more detailed discussion of the ex- 
plicit relation between effective masses and entrainment 
in a two-fluid model we refer the reader to [34]. In the 
context of superfluid mixtures the importance of the in- 
teraction and the entrainment effect has first been recog- 
nized by Andreev&Bashkin [35], although expressed in 
the conceptually more confused orthodox framework of 
superfluidity. Substituting (65) together with (68) and 
(67) into (20), we can now relate the "generalized pres- 
sure" VP directly to the energy function £ , namely 

£ + ^ = Y, n x^ , (69) 

and with (66) the total differential of ^(^f^Aj-y) is 
found as 

= Y,nxdf -W* XY dA XY . (70) 

X,Y 

We can further express the hydrodynamic force density 
(14) more explicitly as 

fl=n x (dt+v x ■ V) pX+nxV^-Y^ 2a XY A XY Vv Xj , 

Y 

(71) 



4 Sometimes also referred to as "drag" in the superfluid literature. 
But as pointed out in [32], this is rather misleading, as entrain- 
ment is a purely non-dissipative effect, whereas "drag" in physics 
usually refers to a resistive drag. 



and for the conserved hydrodynamic energy density (26) 
we find 

£ H = 5>*r^+£-£c^A| y . (72) 

X X,Y 

This relation can be used to clarify the physical meaning 
of the thermodynamic potential £. One might have ex- 
pected to find the total energy density simply as the sum 
of kinetic energies plus £. It is to be noted though that 
E H , which represents the Hamiltonian H H (^x,P fY ) of the 
system, is naturally a function of the fluid momenta pF 
as opposed to the velocities. Similarly it turns out that in 
order to find the actual "internal energy", we have to con- 
struct the thermodynamic potential that depends on the 
relative momenta instead of Axy. We therefore define 
the "entrained" relative momenta J x Y as 

J XY = 2a XY A XY , (73) 

representing the momentum exchange between con- 
stituents X and Y due to entrainment, namely by using 
(67) the momentum density of the constituent X can be 
written as 

nxP* = n x m x v x - ^ J* Y ■ (74) 

Y 

Using this definition of J x Y , the first law (66) now takes 
the form 

d£ = Y,^ dnx + \Y, J XY dA XY , (75) 

X,Y 

We can therefore introduce the internal energy density 
£ as the Legendre transformed (with respect to the mo- 
menta J XY ) of the energy function £, namely 

£(n x ,J XY )^£-l-J2 jXY ■ AxY > ( 76 ) 

X.Y 

with the associated total differential 

d£ = j2^ dn x-lJ2 AxYdjXY ■ ( 7? ) 

X,Y 

We note that £ and £ only differ in systems where the 
entrainment effect is present. Traditionally the quantity 
£ is what one might call the actual "internal energy" den- 
sity, which is a function of the momenta, while the con- 
jugate thermodynamic potential £ does not seem to have 
a well established name in the literature. We see that in 
terms of the internal energy £, the total energy density 
(72) does indeed have the expected form of "kinetic plus 
internal" energy, namely 

En=J2 mXnx Y +i " (?8) 

X 
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C. Entropy and temperature 

As noted earlier, entropy can be included quite nat- 
urally in this framework as a constituent. The corre- 
sponding density and current are n s = s and n s = sv s in 
terms of the entropy density s and its transport velocity 
v s . The entropy is naturally mass- less, i.e. m s = 0. The 
thermodynamically conjugate variable to the entropy (its 
"chemical potential") is the temperature, i.e. fi s — T, so 
(66) can be written as 

d£ = Tds + ^2^ dn x + ^^a XY dA XY . (79) 

X^s X,Y 

The thermal momenta p$ = G>o and p s — of the en- 
tropy constituent are found from (67) and (68), namely 

© = -E^T~ A ^< ( 8 °) 

Y 

-9 = T + v s -@. (81) 

We see that although the entropy has zero rest mass, it 
can acquire a non-zero dynamical momentum due to 
entrainment. This can also be interpreted as the entropy 
having a non-zero "effective mass". The hydrodynamic 
entropy force density /| and energy change rate g s de- 
fined in (14) and (16) yield 

/ S H = sVT+s (8 t + iv V) 0-^2a sY A^V % - ,(82) 

g s =v s -fl + (T + v s -®)r s . (83) 

We see that the temperature gradient is a driving force 
of the entropy constituent, as would be expected. We 
also recognize the term TT S in the expression of the en- 
ergy transfer rate 5 s , which represents the heat creation 
"TdS". 

VI. APPLICATIONS 

A. Single perfect fluids 

As the first application of the foregoing formalism, we 
consider a single perfect fluid consisting of several co- 
moving constituents. This multi-constituent fluid is de- 
scribed by the densities nx which move with a single 
velocity vx = v, and so the currents are nx = nx v. Ob- 
viously all the relative velocities vanish in this case, i.e. 
A^y = 0, and therefore there is no entrainment. Here 
we will explicitly write the entropy with its density s, and 
we do not include it in the constituent index set labelled 
byX, i.e. I/s. The Lagrangian (65) for this system is 

A h = 2Zm x n x — -£ (s,n x ) , (84) 

and the energy and pressure differentials (66) and (70) 
simply read as 

dS = Tds + ^2i^ dn x , and dP = sdT + d^ , 

(85) 



where in the case of a single fluid, the generalized pressure 
\t simply reduces to the usual fluid pressure P. The fluid 
momenta (67) and (68) are 

v 2 

p x =m x v, and -p% =p? +m x —, (86) 

while for the entropy constituent we have with (80) and 
(81): 

= 0, and - O = T . (87) 

The explicit expression for the force densities (11) and 
energy transfer rates (16) are found as 

J* = n x m x (d t + v- V) v + n x V ^ + m x T x v , (88) 

g x = v f +r x f -m xV -r x , (89) 
f s =sVT, (90) 

5 S - v ■ f s + rr s , (91) 

If we allow for an external force / ex t and energy exchange 
rate g C xt, the equations of motion (18) of the system are 

/ s + Yl f x = ' and -9 s + 2 f = - 9 ° xt • ( 92 ) 

Inserting (88)-(91) and using mass conservation (57), we 
find the explicit equations of motion 

(d t +v- V)v + -VP = -/ ext , (93) 
P P 

TTs + J^^rx = ffext - V ■ /ext , (94) 

where we have used the thermodynamic relation (85) in 
order to rewrite the momentum equation in the famil- 
iar Euler form. The energy equation expresses the heat 
creation TT S by chemical reactions Jx- For an isolated 
system, where / ex t = and g cx t — 0, that entropy can 
only increase due to the second law of thermodynamics, 
so f s > 0. From (94) we therefore obtain a constraint on 
the direction of the chemical reactions, namely 

r x ^ < . (95) 

If we consider for example the case of two constituents 
of equal mass, so that the mass-conservation (57) implies 
A + r-2 = 0, then this constraint now reads as 

AO* 1 -M 2 ) < 0, (96) 

which shows that chemical reactions only proceeds in the 
direction of the lower chemical potential as would be ex- 
pected. 

B. "Potential vorticity" conservation: Ertel's 
theorem 

We now consider the case without chemical reactions, 
in which the general perfect fluid discussed in the fore- 
going section can be described effectively as a fluid con- 
sisting only of a single matter constituent and entropy. 
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In this case we can show that the vorticity is generally 
not conserved, but that a weaker form of the vorticity 
conservation still holds. The fluid is described by the 
particle number density n, the mass per particle m and 
a comoving entropy density s. Mass conservation (57) in 
this case reduces to f = 0. If we assume the system to be 
isolated, i.e. / + / s = 0, then the only force per parti- 
cle (33) acting on the matter constituent is the "thermal 
force" (90), namely 

-/ H = -sVT, (97) 
n 

where s = s/n is the specific entropy. If s is constant 
everywhere, then this "thermal force" is conservative, i.e. 
5 = and by (41) the circulation is therefore conserved. 
In the non-uniform case, however, we find 

/ l=-i s dT > (98) 

at JOT, JdT, 

which vanishes only if we integrate along a path <9£ that 
lies completely in a surface of constant s. We can also see 
this in the Helmholtz formulation, namely by applying an 
exterior derivative to (97), one obtains 

d$=-ds/\dT, i.e. Vx#=-VsxVT, (99) 

and it follows therefore from (36) that the vorticity is 
no longer generally conserved in this case. However, 
the quantity ds A d$, or its equivalent dual expression 
Vs • (V x 5), still vanishes identically. Based on this ob- 
servation we construct the "potential vorticity" 3-form Z_ 
as 

Z = dsAw, (100) 
and the dual scalar z is 

Zijk = ze ijk , and z = ^e ljk Z ijk = Vs- (V x p) , 

(101) 

where the last expression was found using (30) . The evo- 
lution of the potential vorticity 3-form Z_ under transport 
by the fluid is 

(d t + £ v )Z = d[(dt + £ v )s\ Aw, (102) 

and therefore Z_ is conserved for isentropic flow, i.e. if 

T s = (d t + £ v )s = Q. (103) 

The dual version of (102), namely the conservation of the 
scalar z is then found as 

d t z + V • (zv) = . (104) 

Traditionally this conservation law is often expressed in 
terms of the scalar a = z/p, which then results in the 
following form of the conservation law: 

{dt + v- V) a = 0, (105) 

which is generally known as "ErteFs theorem" [36, 37]. 



C. Thermally conducting fluids 

We have so far only considered perfect fluids, which are 
perfect heat insulators as the entropy is strictly carried 
along by fluid elements and no heat is exchanged between 
fluid elements. It is quite straightforward to extend this 
to thermally conducting fluids simply by dropping the 
assumption that the entropy flux is bound to the matter 
fluid flow, i.e. we just have to allow v s ^ v, where v s and 
v are the velocities of the entropy fluid and the matter 
fluid respectively. For simplicity we consider only a sin- 
gle matter constituent, described by its particle number 
density n, which by (57) is automatically conserved, i.e. 
r = 0. 

From the general expressions (81) and (80) we see that 
the "entropy fluid" acquires a non-zero momentum due 
to the interaction with the matter fluid, via entrain- 
ment. However, this aspect does not usually seem to be 
taken into account in the traditional description of heat- 
conducting fluids (e.g. see [23]). The aim of the present 
section is only to show how to recover the standard equa- 
tions for a heat-conducting fluid, and we therefore simply 
assume the entrainment to be negligible, i.e. a = 0. It 
is certainly an interesting question if this neglect of en- 
trainment is physically justified in all cases. With this 
assumption, the force density (82) and energy rate (83) 
of the entropy reduce to 

/ s = sVT , and g s = v s ■ f + T T s . (106) 

As in the (isolated) perfect fluid case discussed previ- 
ously, the equations of motion are again / s + / = and 
g s + g = 0. This time, however, one force density, / s say, 
can be specified by the model due to the increased num- 
ber of degrees of freedom, so we set it to / s = fa, where 
fit is a resistive force acting against the entropy flow. We 
obtain the Euler equation in the same form as in (93), 
but now the energy equation takes the form 

Tr s = (v- v s ) ■ f R . (107) 

By the second law of thermodynamics, namely r s > 0, 
we can constrain the form of the resistive force fa to 

fR = -V (v s - v) , with r] > , (108) 

i.e. the friction force acting on the entropy fluid is always 
opposed to its flow relative to the matter fluid. Obviously 
the value of the resistivity rj is not restricted to be a 
constant but will generally depend on the state of the 
system. Following the traditional description (e.g. [23]) 
we introduce the heat flux density q relative to the matter 
fluid as 

q = Ts(v s - v) . (109) 

By combining this with (106) and (108), we see that the 
heat flux current is constrained by the second law to be 
of the form 

Ts 2 

q = -KVT, with k= >0, (110) 

n 



12 



where k is the thermal conductivity. With (109) we can 
express the velocity of the entropy fluid v s in terms of 
the heat flux q and the matter velocity v, so the entropy 
creation rate r s can be expressed as 



r s = d t s + v • («» + |) . 



(in) 



We further find for the hydrodynamic energy flux vector 
Q H of (26): 



Qn = 22(-Po)nx = (fi + m—)nv + sTv s 



nv [ m— + u + sT ) + q , 



(112) 



where the last equality was found using (109). We in- 
troduce the specific enthalpy as w = ji + sT, and using 
the first law 5 , namely dP — ndp + s dT, we find the total 
variation of the specific enthalpy as 



dw = Tds + -dP . 

n 



(113) 



and so we recover the standard expression (e.g. cf. [23]) 
for the energy flux: 



(114) 



Q H = nv ( m— + w ) + q . 



D. The two-fluid model for superfluid 4 He 

We now consider the example of superfluid 4 He at a 
non-zero temperature T. Let n be the number density of 
4 He atoms and s be the entropy density. The 4 He atoms 
move with a velocity v, while the entropy (carried by a 
thermal gas of excitations such as phonons and rotons) 
transports heat without friction (i.e. fu = 0) at the 
velocity »n, so the relative velocity is A = dn — v. In 
this context the entropy fluid is often referred to as the 
"normal fluid" as opposed to the superfluid mass flow. 
The two transport currents, namely that of 4 He atoms 
and of entropy, are respectively 



n = nv , and s = sdn. 



(115) 



The 4 He atoms have mass m, so the mass density is 
p = nm, and the hydrodynamic Lagrangian density (65) 
reads as 



A H = -nmv 2 — £(n, s, A 2 ) . 



(116) 



6 In the absence of entrainment the entropy fluid does not carry 
momentum, therefore the matter fluid defines a unique frame in 
which the stress tensor (23) is purely isotropic. In this case the 
generalized pressure ^ is identical with the usual perfect fluid 
notion of the pressure P. 



where the energy function £ determines the first law (66) 
as 



d£ = udn + Tds + adA 2 



(117) 



which defines the chemical potential n of 4 He atoms, the 
temperature T and the entrainment a. The conjugate 
momenta (67), (68) of the 4 He atoms are 



2a 

p = mv H A , 

n 

1 2 

-po = \i — -mv +v-p, 



(118) 
(119) 



while for the entropy fluid Eqs. (80) and (81) yield 



= -^A. 

s 



(120) 

-Go = T + v N -&. (121) 

The conservation of mass (57) implies 

r = d t n + X7 ■ n = 0. (122) 

In the absence of vortices, there are no direct forces acting 
between the two fluids, so the equations of motion in the 
absence of external forces (i.e. / ex t = 0) are simply 



/ = / H = and f = . 



(123) 



The energy equations are g = and <? N = g ex t, and with 
(83) this leads to 

-Sext - r s (9 + «n • ©) = -Tr s , (124) 

where we have inserted (121). We see that this equation 
describes the rate of entropy creation by an external heat 
source, namely 



d t s + V • (sv N ) = ^ cx t- 



(125) 



As discussed in Sect. IV B, the superfluid 4 He is (locally) 
irrotational, i.e. 

wn = 2V [i p j] =0,^W = Vxp = 0. (126) 

Using (14), the equation of motion (123) for the super- 
fluid therefore reduces to 

dtp - Vp = , (127) 
and with the explicit momenta (119) and (118) this yields 

d t (v + eA) + V (jl +^v 2 +ev- A^j = 0, (128) 

where we introduced the entrainment number e and the 
specific chemical potential ju as 



2a „ a 

e = — , and n = — 
p m 



(129) 
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The entropy fluid is governed by the momentum equation 
/ N = 0, and with (82) and the entropy momenta (121) 
and (120), we find 

(d t +v N -V) (^*) - VT + y Aj-Vvj, + ^r s A = . 

(130) 

The two equations (128) and (130) represent the "canon- 
ical" formulation of the two-fluid model for superfluid 
4 He. These equations do not seem to bear any obvious 
relation to the "orthodox" formulation of Landau's two- 
fluid model found in all textbooks on the subject (e.g. 
see [23, 24, 31]). Nevertheless, these equations are equiv- 
alent to the orthodox framework, as we will show now, 
but it is important to note that the orthodox formulation 
is based on a rather unfortunate confusion between the 
velocity and momentum of the superfluid which is inher- 
ent in the historic definition of the "superfluid velocity" 
by Landau. 

We demonstrate the equivalence of these formulations 
by explicitly translating the canonical formulation into 
the orthodox language. The starting point of Landau's 
model is the statement that the "superfluid velocity" is 
irrotational. We write vs for the "superfluid velocity", 
which is not to be confused with the actual velocity v of 
4 He atoms, so the starting point is 



Vxi/ s = 0. 



(131) 



From the general discussion about vorticity conserva- 
tion in Sect. IV A and its particular role in superfluids 
(Sect. IV B) we have already seen that contrary to the 
momentum vorticity W = V x p, the velocity-rotation 
V x v is generally not conserved by the fluid flow, and 
in particular not in the presence of more than one fluid as 
is the case in superfluid He at T > 0. The only possible 
interpretation we can give ug in order for the constraint 
(131) to be consistent with hydrodynamics and to remain 
true for all times is that it is really the rescaled super- 
fluid momentum p, so the "key" to our translation is the 
ansatz 



P_ 

m 



(132) 



While this would be equivalent to the fluid velocity in a 
single perfect fluid, as seen in (86), this has no interpre- 
tation as the velocity of either the mass or the entropy 
in the case of the present two-fluid model as we can see 
in (118). Therefore we call vs a pseudo velocity, as it is 
a dynamic combination of both fluid velocities governed 
by the entrainment a between the superfluid 4 He and 
its excitations. With the explicit entrainment relation 
(118) we can now express the velocity v of the 4 He fluid 
in terms of the pseudo-velocity vs and the normal-fluid 
velocity «n as 



v = (1 -e) 1 (i/ s - e^n) : 



(133) 



where we used the definition (129) of the entrainment 
number s. With this substitution, the total mass current 



p, which is equal to the total momentum density J 

H aS 

seen in (63), can be written in the form 



J H = pv 



1 



-ep 



1 



(134) 



which suggests to introduce a "superfluid density" gs and 
a "normal density" as 



Qs 



P A _ -£P 

, and 0n 



1-e 



1-e 



(135) 



such that total mass density p and mass current p = J H 
can now be written as 



P = gs + £>n , and J H = £>s^S + £»n«n • 



(136) 



It is now obvious that this split is completely artificial, 
and £?n and gs are only pseudo densities, as they do not 
represent the density of any (conserved) physical quan- 
tity and are not even necessarily positive. In fact neither 
of the two pseudo-densities and currents are conserved 
individually, contrary to the physical currents (115). We 
note that even Landau warned against taking too liter- 
ally the interpretation of superfluid 4 He as a "mixture" of 
these two (pseudo-) "fluids" [23]. Contrary to the artifi- 
cial orthodox split, however, the separation into entropy 
fluid and 4 He mass flow is physically perfectly meaning- 
ful, and the superfluid can be regarded as a two-fluid 
system in the literal sense in the canonical framework. 
The pseudo "mass density" £>n, which the normal fluid 
seems to carry in the orthodox description is due to the 
fact that entrainment provides the entropy fluid with a 
non-vanishing momentum (120) in the presence of rela- 
tive motion, even though it does not transport any mass. 
This lack of careful distinction between mass current and 
momentum leads to the paradoxical picture of the "super- 
fluid counterflow": for example, in the simple case of heat 
flow through a static superfluid, the normal fluid associ- 
ated with the heat flow carries a pseudo mass-current 
Pn^n- But because there is no net mass current there 
has to be some superfluid "counterflow" of pseudo mass 
current gs^s — — Pn^n- This apparently strange behav- 
ior is solely due to an awkward choice of variables and a 
loss of direct contact between the quantities used in the 
orthodox description and the actual conserved physical 
quantities of 4 He. 

Further following the traditional orthodox framework, 
we define the relative (pseudo-)velocity w as 



W = t)N — t'S : 



(137) 



which, using (133), can be expressed in terms of A as 

w = (1-e) A. (138) 

In order to relate the canonical thermodynamic quan- 
tities to the orthodox language, we follow Khalatnikov 
[31] and Landau [23] and consider the energy density in 

the "superfluid frame" Kq, which is defined by i/g = 0. 
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In this frame, the momentum density ' expressed in 
(136) is 

4 0) = PN v { ° ] = g N w = -2a A , (139) 

and the transport velocity v of the superfluid 4 He atoms 
in this frame can be expressed using (140) as 

v W=v-v s = ^vW = ±jM. (140) 

The hydrodynamic energy density E H of the fluid system 
is given by (72), which reads in this case 

E H = \pv 2 +£ -2a A 2 , (141) 

and using the previous translations together with the first 
law (117), we can write the total variation dE^ of the 
energy density in Kq as 

dE ( ° } =Tds + Jlsdp + w dJ { ^ , (142) 
which defines the "superfluid chemical potential" Jls as 

£s =£- \{v-u s f. (143) 

Using these quantities, the canonical equation of motion 
(128) can now be translated into the orthodox form as 

^s + V^+ms)=0. (144) 

One can equally verify that the generalized pressure, de- 
fined in (69), is expressible in terms of the orthodox quan- 
tities as 

* = -£ + pJl+sT= -E^+Ts + pjis+w-jW , (145) 

in exact agreement with the expressions found in [23, 
31]. For the remaining momentum equation, the total 
momentum conservation (22) is traditionally preferred 
over the equation of motion (130) of the entropy fluid. 
We therefore conclude this section by the appropriate 
translation of the stress tensor (23) into the orthodox 
language. The canonical expression for the stress tensor 
of superfluid 4 He is 

Tjj 7 ' = n l p> + s i e J + * g 13 , (146) 

and inserting the previous expressions for the explicit mo- 
menta and the translations to orthodox variables, one can 
write this in the form 

T% = ?s 44 + £N 4< + * 9 ij , (147) 

which concludes our proof of equivalence between canon- 
ical and orthodox description. 



E. A two-fluid model for the neutron star core 

Here we consider a (simplified) model for the matter 
inside a neutron star core, which mainly consists of a 
(charge neutral) plasma of neutrons, protons and elec- 
trons. We focus on superfluid models, in which the neu- 
trons are assumed to be superfluid, which allows them 
to freely traverse the fluid of charged components due to 
the absence of viscosity. As discussed in Sect. IV B, this 
also implies some extra complications due to the quanti- 
zation of vorticity into microscopic vortices. Here we are 
interested in a macroscopic description, i.e. we consider 
fluid elements that are small compared to the dimensions 
of the total system, but which contain a large number of 
vortices. On this scale we can work with a smooth aver- 
aged vorticity instead of having to worry about individual 
vortices. One effect of the presence of the vortices will be 
a slight anisotropy in the resulting smooth averaged fluid 
[22, 38, 39], which can be ascribed to the tension of vor- 
tices, and which we will neglect here for simplicity. The 
second effect of the vortex lattice is that it allows a direct 
force between the superfluid and the normal fluid, medi- 
ated by the respective vortex interactions, and which is 
naturally described in the context of the two-fluid model 
as a mutual force. The model assumptions used here 
are fairly common to most current studies of superfluid 
neutrons stars (e.g. see [34, 40-42]). 

The model therefore consists of comoving constituents 
X e {e, p, s}, corresponding to the electrons, protons and 
entropy, and we will label this fluid with V. The second 
fluid consists only of the superfluid neutrons, i.e. X = n. 
Charge conservation implies 

r c = r p , (148) 

and for simplicity we will assume local charge neutrality, 
i.e. 

n c = n p . (149) 

We assume the electrons and protons to be strictly mov- 
ing together in this model (i.e. we consider timescales 
longer than the plasma oscillation timescale), so we can 
neglect electromagnetic interactions altogether. Another 
physical constraint is baryon conservation, i.e. we must 
have 

r n + r p = o , (150) 

and together with mass conservation (57), this leads to 
the requirement 6 

m n = m p + m c = m. (151) 



This relation is of course not exactly satisfied in reality, which 
shows a well-known shortcoming of Newtonian physics: mass has 
to be conserved separately from energy. 
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We can therefore write the mass densities of the two fluids 

as 

p n ~mn D , and p c = mn p . (152) 

The first law (66) of this model reads as 

d£ = Tds + p n dn n + p c dn c + pP dn p + a cn dAj n 

+a pn dA 2 pD + a sn dA 2 sn . (153) 

Obviously there is only one independent relative velocity 
A, namely 

A = v c - v n = A on = A pn = A sn , (154) 
and we define the total entrainment a as 

a = a cn + a pn + a sn . (155) 

In the case of the neutron star model, we are obviously 
also interested to include the effects of gravitation. We 
can therefore not assume the system to be isolated and 
we include the effect of the gravitational potential $ as 
an external force. The minimal equations of motion (18) 
therefore read as 

P + f c = -pV$ , and <? n + g c = p • V$ , (156) 

where the force and energy rate of the 'c'-fluid are natu- 
rally given by f c = fP + f c + / s and g c = g? + g c + g s . 
With (148) and (150) we can write the respective force 
densities more explicitly as 



r = + (157) 

f c = r H -r n (p c + P v) + r s ®, (158) 

where we naturally defined /£ = /h + /£ + /h • Simi- 
larly we can write the energy rates (16) as 

g n = v D ■ r H - r nP » , (159) 

9 C = v c ■ fi + r n ( P e + pg) - r s e . (ieo) 



Because the gravitational acceleration is the same for all 
bodies (i.e. fluids), we can now simply absorb the ef- 
fect of the gravitational potential into the definition of 
"extended" forces / and energy rates g which simply in- 
corporate the respective gravitational force density and 
work rate, i.e. we define 

f = f X +PxV$, (161) 

fl = /^+PxV$, (162) 
<f = gX+pxVx-V®. (163) 

With these redefinitions, the minimal equations of mo- 
tion (156) again take the form of an isolated system, i.e. 

/ n + / C = 0, and ^+^ = 0, (164) 



while for (157)-(160) we obtain exactly the same form, 
just for all forces and energy rates replaced by their "ex- 
tended" version. Using the foregoing equations, we ob- 
tain 

?h = + r n p c - r s & , (165) 

and therefore 

F = -v c ■ ?° - r n [v c ■ ( P n - P c ) - pi] - r s e . (166) 

Substituting this and the "extended" version of (159) into 
the energy-rate equation (164), we find 

rr s = A-T H + r n lrt-p c -p p + v c -(p»-p c - P v)} , 

(167) 

where we have used the explicit form (81) of 6o- In addi- 
tion to the external force, the two-fluid model allows one 
to prescribe one of the fluid force densities. In the present 
case it is most convenient to specify the "extended" hy- 
drodynamic force / H on the neutrons. As this force can 
only originate from the second fluid, we will refer to it as 
the mutual force / mu t, so we set 

/h - /mut • (168) 

Substituting the explicit conjugate momenta (67) and 
(68), we obtain the final expression for the entropy cre- 
ation rate (167) as 

TT S = A 

J mut + r n /3 . (169) 

The first term on the right hand side is the work done by 
the mutual force, and the second term is the entropy cre- 
ated by beta reactions between the two fluids, for which 
the term "transfusion" has been coined [32]. The devi- 
ation from beta equilibrium characterized by (3 is found 
as 

p^pV + f-^- l - m (l-^ A 2 , (170) 

where the last term gives the correction to the chemical 
equilibrium due to relative motion A of the two fluids. 
The second law of thermodynamics for an isolated system 
states that entropy can only increase, i.e. r s > 0. In 
order for this to be identically true in (169), the mutual 
force /mut and the reaction rate T n have to be of the 
form 

T n = Ef3, with 5 > 0, , 171 - 

/ mut = t)A|cxA, with 1] > , 1 > 

where k is an arbitrary vector characterizing a non- 
dissipative Magnus-type force orthogonal to the relative 
velocity. Further substituting the conjugate momenta in 
the expression for the hydrodynamic force densities (14), 
we find their explicit form 

fl=n n (d t + v n -V) (mv n + ^Aj +n n \7p n + 2aA 3 V< , 

(172) 
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f c H =n p (d t + v c -V)[mv c ^ ^A)+n p V( M p +^ c ) 



/2a sn 

-2aAjVvt - s(d t + v c ■ V) A ) + sVT. (173) 



V s 

We now make the simplifying assumption that we can 
neglect the entrainment of entropy, i.e. we assume that 
all the entrainment between the two fluids is due to the 
neutron-proton and neutron-electron contributions, so we 
set a sn = 0, which implies © = 0. Using (67) we find 

p c +p p -p n = m(l-e n -e c ) A, (174) 

where we have defined the entrainment numbers 



_ 2a , _ 2a 
E n = — , and e c = — 

Pn Pc 



(175) 



Putting all the pieces together, we obtain the momentum 
equations (168) and (165) in the form 

(d t + v n • V) (v n + e n A) + V (Jl n + *)+e n A 3 V< = — / mut , 

Pn 

(176) 

{d t + v c -V)("y c -£ c A) + V(M c + $)-£cA J V^ + -VT 

Pc 

= -— /mut + (1 - £c - £n) m— A. (177) 

Pc Pc 

with the specific chemical potentials p, n = /i n /m and 
Jl c = {pP + p c )/m. 
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APPENDIX A: EVALUATION OF CONVECTIVE 
VARIATIONS 



We write the particle flowlines as 

x 1 = x l (a, t) , 



(Al) 



where the "particle coordinates" a 1 are used to label in- 
dividual particles and can be taken, for example, to be 
their initial position, i.e. 



This introduces a time-dependent map (or "pull-back") 
between the "material space" a 1 and physical space x l , 
and the associated Jacobian matrix J is 



V = 



dx l 



(A3) 



We consider the variations of fluid variables induced 
by active infinitesimal spatial displacements £ l (x,t) and 
temporal shifts t(x, t) of the fluid particle flowlines (Al), 
namely 

x H (a,t') =x i (a,t)+C(x,t), and t' = t + t(x, t) . 

(A4) 

We note that the transformation (A4) not only shifts 
flowlines in space, but also in time. A physical quantity 
of the flow, Q(x,t) say, is changed to Q'(x',t'), and we 
define the corresponding Eulerian and Lagrangian varia- 
tions as 7 

5Q= Q'{x,t)-Q(x,t), (A5) 
AQ= Q'(a,t')-Q(a,t) = Q'(x',tf)-Q(x,t). (A6) 

By expanding AQ to first order using the definition (A4) 
of x l and t' we find the relation 



AQ = SQ + e VjQ(x, t) + r d t Q(x, t) . 



(A7) 



Let us consider the induced (first order) variation of the 
velocity v % = dtx % (a,t), namely 



v H {a,t') = d t ,x' i {a,t') = d t ,x i {a,t) + d t tt{a,t) 
dt 

= d t x*(a,t) — +d t C(a,t) 

= v l {a,t) - v l d t r(a,t) + d t ?(a,t) , (A8) 



which by (A6) corresponds to the following Lagrangian 
variation of the velocity: 

Av l = [dtC + i/Vzf ] - [v l d t r + V/r] , (A9) 

and with (A7) the Eulerian variation is found as 

Sv* = [d t C + v l ViC - ~ [dt (v*t) + v*v l V lT ] . 

(A10) 

From the conservation of mass one can derive an expres- 
sion for the particle density n in terms of the Jacobian 
(A3), namely 



n(x, t) = 



n {a) 
detj 



(All) 



a 1 = x l (a,0). 



(A2) 



7 Contrary to the Eulerian variation, the Lagrangian variation can 
be defined in different (non-equivalent) ways. The definition used 
here is based on comparing the quantity Q in different points 
by parallel-transport. Another common definition (e.g. see [11, 
20]) consists in using the Lie-transported quantity instead. Both 
definitions are equivalent for scalars but differ for vectors and 
higher order tensors. 
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where no (a) = n(a,0) is the initial density at t = 0. 
Using (A3) , the change of the Jacobian matrix J induced 
by the flowline variation (A4) can be found as 



dx H (a,t') _ dx'fat) 



dx t {a 1 t) dx l (a,t) dt 
da? 



+ 



da? 



da? 



dt 



da? 



with the resulting Lagrangian variation (A6) expressible 
as 



(A13) 



The derivative of a determinant det A with respect to a 
matrix element Aij is given by 



ddctA 
dAa 



= det(A) (A 



-1\ % 3 



(A14) 



and therefore we can write the Lagrangian variation of 
the Jacobian determinant as 



A (det J) = det(J) (J' 1 )'. AJ l j 



(A15) 



The flowline variation (A4) therefore induces the La- 
grangian change of the Jacobian 

A (det J) . , , A N 

-L rf i = V l€ '-«'V I r. (A16) 

Using (All), the induced density variation is therefore 
found as 



An = -nVi? +nv l Vir, 



(A17) 



and with (A7) the corresponding Eulerian expression is 
found as 



Sn 



-V; (<') + [nv l X7 lT -Td t n] 



(A18) 



By combining the results for velocity and density varia- 
tions we find the variations of the current n % = nv l as 

An l = [ndtC(x,t) + n l Wie - n 4 V^] - n l 8 t T , (A19) 
«5n 4 =[na t f( a; ,t) + n z V ; f-V i (n^)]-5 t (nV) . (A20) 



APPENDIX B: NOETHER IDENTITIES OF THE 
VARIATIONAL PRINCIPLE 

In addition to the flowline variations considered so far, 
we will now also allow for metric variations Sgij. Al- 
though we only consider Newtonian physics here, there is 
a-priori no reason to restrict ourselves to flat space. Most 
importantly, however, including metric variations allows 



us to obtain the form of the stress tensor T„ and the as- 
sociated momentum conservation (22) directly from the 
variational principle as a Noether identity, as opposed to 
constructing it from the equations of motion as we have 
done in Sec. III. Therefore we extend the variation (4) 
of the Lagrangian to 



dgij 



(Bl) 



(A12) Next consider the density change Sn x induced by a metric 



variation 



at constant flowlines, i.e. constant J' 



First we note that we can express the Jacobian as 
detJ=e ijk J i 1 J j 2 J k 3, 



(B2) 



and using (A14) the variation of the volume form 
e ijk = \fg [ijk] induced by metric changes is expressible 
as 



Seijk 

Therefore we have 

ddctj 



dgi 



■det(JV 



(B3) 



(B4) 



3 



and using (All) and (A18) we can write the variation 
of the density induced by spatial displacements £ and 
metric variations Sgij as 



Sn = -V; (<') - ^ng' 13 Sg k 
An = -nX7,£. 1 - ^ng 11 5g l0 , 



(B5) 
(B6) 



where we have used the fact that with our definition of 
the Lagrangian variation (A7) we have 

&9ij = $9ij + ^l9ij = s 9ij , (B7) 

as the metric is by definition constant under parallel 
transport. A metric change with fixed flowlines does not 
change the local velocity v % , therefore the current varia- 
tion can be written using (B5) and (A20) as 

Sn* = [ndtC {x,t) + n l V£-V l („*£')] -I n V j ''Sgij , 

(B8) 

An 1 = [nd t e(x,t)+n l V l e-n*V l C l ] - \n l g l 'S 9l] . 

(B9) 

When allowing for metric variations it is convenient (e.g. 
see [39]) to introduce the "diamond variation" 0A H as 

0A H = -j=S (V? A H ) = 5A H + \K , (BIO) 

such that the variation of the action (3) can now be writ- 
ten as (noting that dV — y/gd 3 x): 



SI 



J OK 



dVdt. 



(BH) 
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Substituting (Bl), (B5) and (B8) and integrating by 
parts, 0A H can be cast in the form 

0A H = - ]T f x Qc + l -T% Sg tJ + V t R l + 8 t R , (B12) 

where the canonical forces fx have the explicit expres- 
sion (11) and we defined the tensor as 



d\ 

ogi 3 



(B13) 



using our earlier definition (20) of the generalized pres- 
sure VP. 

Now consider a common displacement £ of the whole 
system including the background metric, which induces 
a metric change 



Sgij = -2V (i £j) 



(B14) 



where (ij) indicates symmetric averaging, i.e. 
2A(^B^ — AiBj + AjB{. The corresponding Lagrangian 
variations (B9) and (B6) are found as 

An x = 0, (B15) 
A4 - nx{d t e + v x ViC) ■ (B16) 

Substituting this into (Bl), the induced AA H is 

AA H = (j2 ^ 3 - 2^f) + Jl m , (B17) 

where we have used the definition (23) of the momentum 
density J H . It is well known that contrary to the fully 
covariant Lagrangian for relativistic hydrodynamics (e.g. 
[13]), the Newtonian Lagrangian is not strictly Galilean 
invariant under boosts. This is due to the velocity depen- 
dence of the kinetic energy, as can be seen in the explicit 
form (52). We can therefore only demand strict invari- 
ance, i.e. AA H = 0, for time-independent displacements, 
namely d t £ = 0, which leads to the Noether identity 



dg 



(B18) 



The left-hand side is manifestly symmetric in i and j, 
therefore we see that 



(B19) 



and we can now write the (symmetric) stress tensor (B13) 
explicitly as 



(B20) 



This tensor is identical to the expression (23) found ear- 
lier by construction from the equations of motion. It 
remains to be shown however, how the momentum con- 
servation law (22) is directly obtainable Noether 
identity from the variational principle. Using (B17), (A7) 



and (B12) we can explicitly express the diamond varia- 
tion as 

0A H = -(d t J 3 ) - V ; (A H £') + d t (4 £) , (B21) 

which has to be identical to the expression (B12) for a 
common displacement £ of the whole system, which after 
some partial integrations takes the form 

<> A h = (- E f Xj + V < T h ) + V,(...)' + dt(-) ■ 

(B22) 

The requirement that the previous two expressions have 
to be identical (up to divergences and time derivatives) 
leads to the Noether identity 

dtJi + VjTH =f eiA , (B23) 

which is the momentum conservation law (22). 

APPENDIX C: GALILEAN INVARIANCE OF £ 

In this section we show that requiring the conjugate 
momenta p§ and of (53) to transform as (55) under 
Galilean boosts (54) implies that the internal energy £ 
has to be Galilean invariant. We assume that £(nx,nx) 
transforms into £'(nx,n x ) under a Galilean boost, where 



n' x = n x +n x V . 



(Cl) 



Therefore the conjugate momenta (53) in the frame mov- 
ing with speed —V are of the form 

-PS'= \ m x v x +rr?v x .V+\m x V* + ^-, (C2) 



m x vx 



m V 



d£' 
dn' x ' 



(C3) 



Using (53) to eliminate all terms containing vx , we arrive 
at 



-pl'=-^ + V-^ + \m x V 2 + 



d£' d£ Tr d£ 
+ V 

dn x n x n x 



jp'=p + m x V + 



d£ d£' 



dn 



x 



dn' x 



(C4) 
(C5) 



By comparing with the required transformation proper- 
ties (55) we see that a necessary and sufficient condition 
for this is the vanishing of the terms in brackets in (C4) 
and (C5). We can rewrite the partial derivatives of the 
energy function as follows 



8£' d£' dn x 



dn x dn x dn' x 



and 
d£' 



dn 



x 



d£' 



dn 



•x 



d£' dn x 



dnx dnx 



d£' 
dn x 



d£' 



dn 



•x 



V 



(C6) 



d£' 
dn x ' 
(C7) 
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Inserting these identities into (C4) and (C5), the invari- 
ance requirement can be expressed as 



d£ 



dnx 



d£' 



dn x 



and 



d£ 



d£' 



dnx dnx 



for all A". 



(C8) 

therefore £' can only differ from £ by a constant, which 
is unimportant because the absolute value of the energy 
scale is arbitrary. This shows that energy function £ has 
to be Galilean invariant under the above assumptions. 



APPENDIX D: NEWTONIAN LIMIT OF THE 
RELATIVISTIC LAGRANGIAN 

As shown in the relativistically covariant framework 
by Carter [13], the equations of motion for conducting 
multi-constituent fluids can be derived from a covariant 
Lagrangian density of the form 



Acov = -PC 2 



(Dl) 



where the scalar p is now the total mass-energy den- 
sity of the system. For simplicity we consider here a 
two-fluid system, as generalizations to more fluids are 
straightforward while making the notation more cumber- 
some. The two fluids, A and B say, are described by the 
two 4-current densities n A , n^, and therefore the scalar 
A CO v(«A' n B) can on ly depend on the three independent 
scalar combinations of these two currents, for example 



n A = ^\/-.</,^"V>'\ ■ 



1 



n B = 



9fJ.v n B n B 



and 



9fJ,f n A n B ■ 



(D2) 



and so generally A cov = A cov (nA, n B ,x). Instead of x we 
can equivalently choose as a third independent quantity 
the combination 



A 2 

A 



^ n A n B ^ ' 



(D3) 



We are interested here only in the purely hydrodynamic 
content of this framework, so we assume a flat space-time, 
i.e. a metric of the form 



ds 1 



g^dx^dx" = -(? dt 2 + dx 1 



(D4) 



with the time-coordinate x° = t and so .goo = — c 2 . 
When taking the Newtonian limit as c^> oo, the met- 
ric becomes singular. The reason for this singular limit 
obviously lies in the fact that a locally Lorentzian the- 
ory reduces to a Galilean invariant theory, therefore the 
Lorentz invariance has to be broken in the limit. As 
the non-invertible metric no longer fully determines the 



space-time, we now have to choose 8 a preferred time co- 
ordinate, t say, in which to take the limit and which will 
reduce to the Newtonian absolute time. 

The relation between the scalar rest-frame particle 
densities nx and the densities n x in the preferred-time 
frame can be expressed from (D2) and (D4): 



1 

n x = - 



2 K) 



n 



■x ' 



-n x 



1 



1 / v x 

'2 \ c 



+ 0(c- 4 ) 



(D5) 

where (nx) 1 = n l x is the spatial part of the 4-current n x 
in the preferred time frame, and the relation to the New- 
tonian 3-velocity vx is simply nx — n x v x . We see from 
this equation that if we choose the densities n x to rep- 
resent the Newtonian particle number densities indepen- 
dent of c, then in the limit we find 



lim nx — n x . 



(D6) 



We further note that the quantity A introduced in (D3) 
reduces to the relative velocity in the Newtonian limit, 
namely 

A 2 = (v a -vb) 2 +0((v/c) 2 ) . (D7) 

We now turn to the covariant Lagrangian A cov of (Dl) 
which we can quite generally be written as 

Acov = -{n A m A + n B niB) c 2 - £(n A ,n B , A 2 ) , (D8) 

where the first term represents the rest-mass energy in 
the fluid frame, while £ contains the "equation of state", 
i.e. the internal-energy function of the fluid. When we 
write this in the preferred time-frame using (D5), we ob- 
tain 



-(n° A m A + n^mB)c 2 + \m A n A v 2 A + ]- 



niBn^VQ 



-£(n A ,n B ,A 2 )+O{(v/c) 2 ) 



(D9) 



We see that this Lagrangian obviously diverges in the 
Newtonian limit c — > oo due to the rest-mass energies 
n x mxc 2 . Before we can take this limit, we therefore 
have to renormalize the Lagrangian density by subtract- 
ing a finite counter-term that will make the limit finite. 
The most natural choice is obviously to subtract the 
mass-energy in the preferred-time frame that will deter- 
mine the Newtonian absolute time. We therefore define 
the renormalized Lagrangian density A rcn as 



A ron = A cov + (n° A m A + n B m B ) c 2 



(D10) 



In A rcn we have explicitly broken Lorentz invariance by 
choosing a preferred time frame, and when taking the 
Newtonian limit we obtain the finite Lagrangian 



lim A r 



2n A 



+m B 



2n B 



-£(n A ,n B ,A 2 ), (Dll) 



8 See [15] for a more detailed discussion of this limit and how to 
construct a fully space-time covariant Newtonian framework. 
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which corresponds exactly to the Newtonian hydrody- namic Lagrangian A H of (65). 
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